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\Q • We consider classical Euclidean gravity solutions with a boundary. The boundary 

(N ' 

contains a non-contractible circle. These solutions can be interpreted as computing the 

When the circle is contractible in the bulk, we argue that the entropy of this density matrix 

m ■ 

is given by the area of a minimal surface. This is a generalization of the usual black hole 
entropy formula to euclidean solutions without a Killing vector. 

A particular example of this set up appears in the computation of the entanglement 
entropy of a subregion of a field theory with a gravity dual. In this context, the minimal 
area prescription was proposed by Ryu and Takayanagi. Our arguments explain their 
conjecture. 
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1. Introduction 

Originally the concept of entropy arose from equilibrium thermodynamics. However, 
we know think of entropy as a measure of information. In particular, we can assign an 
entropy to a general density matrix via 

S = -Tr[p\ogp] (1.1) 

By thinking about the thermodynamics of black holes the area formula for gravita- 
tional entropy was discovered [[j],|2],|3|] . Gibbons and Hawking introduced a thermodynamic 
interpretation of euclidean gravity solutions with a U(l) isometry [|J. The idea is that one 
considers Euclidean solutions with prescribed boundary conditions. The boundary condi- 
tions, as well as the solutions, are invariant under a U{1) symmetry!!. These solutions can 
be viewed as describing the computation of the partition function of a quantum theory 
in the classical approximation. In other words, one thinks of the Euclidean gravitational 
action as logZ(/5) = —SE, gr av Then the entropy, obtained as S = —{fidp — l)logZ, is 
equal to the area of the codimension two surface which is a fixed point for the U(l) sym- 
metry in the bulk. Classically, the boundary can be chosen to be any surface where we put 
boundary conditions. It can also be an asymptotic boundary such as the AdS boundary. 





(a) (b) 

Fig. 1: (a) A euclidean solution with a U(l) symmetry is interpreted as computing 
the equilibrium thermodynamic partition function of the gravity theory, (b) We 
consider a euclidean solution with a circle but without af/(l) symmetry. This is 
interpreted as computing Tr[p] for an un- normalized density matrix in the gravity 
theory. This is the density matrix produced by euclidean evolution. 



Here we assume that there is a single U(l) symmetry, otherwise we need to add the corre- 
sponding chemical potentials, etc. 

1 



Interestingly, one can extend the notion of gravitational entropy to situations without 
a U(l) symmetry as follows. 

Let us first consider a general quantum system. Its Euclidean evolution generates an 
un-normalized density matrix 

p = p e J -o (1.2) 

where we considered a general time dependent Euclidean Hamiltonian. We can compute 
the entropy of this density matrix by the "replica trick" . Namely, first notice that Tr [p] 
can be computed by considering euclidean evolution on a circle, identifying Tf = r + 2% 
u. Similarly, we can compute Tr[p n ] by considering time evolution over a circle of n times 
the length of the original one, where the couplings in the theory are strictly periodic under 
shifts of the original circle, H{r + 2-k) = H(t). 
We then can compute the entropy as 

S = - nd n [logZ(n) - nlogZ(l)]\ n=1 = -Tr[p log p] , 

Z(n) = Tr[p«], p^^j-j (L3) 

where now p is a properly normalized density matrix. This involves computing Z(n) and 
then performing an analytic continuation in n. 

Going back to the gravitational context, we can consider metrics which end on a 
boundary. We assume that the boundary has a direction with the topology of the circle. 
The boundary data can depend on the position along this circle but it respects the period- 
icity of the circle. We define the coordinate r ~ r + 2n on the circle. We can then consider 
a spacetime in the interior which is smooth. Its Euclidean action is defined to be log Z(l). 
See fig. 1(b). We can also consider other spacetimes where we take the same boundary 
data but consider a new circle with period t ~ r + 2nn. Their Euclidean action is defined 
to be log Z{n). These computations can be viewed as computing Tr[p n ] for the density 
matrix produced by the Euclidean evolution. See fig. 2. If we are sufficiently diligent, we 
can find these actions, analytically continue in n the corresponding answers and compute 
S as in ( |1.3| ) . This has been explicitly done in []5]U for some examples in three dimensional 
gravity. 

Note that we are implicitly assuming that gravity is holographic. We are imagining 
that setting boundary conditions on some boundary defines the theory and that the interior 



Throughout this paper we set the coordinate length of the initial circle to 2ir. Of course, its 
physical length depends on the metric. 




n = 1 



Boundary 



n = 4 



Fig. 2: Computing the entropy using the replica trick, (a) Euclidean solution 
for n = 1. (b) Solution for n = 4. At the boundary we go around the original 
circle n times before making the identification. We then find a smooth gravity 
solution with these boundary conditions. The curves in the right hand side are 
schematically giving the boundary conditions at infinity. We see that in (b), we 
simply repeat n times the boundary conditions we had in (a). 



geometry is an approximation to the full computation. We do not know how (and whether) 
holography works for general boundaries. Here we only need it to be approximately valid so 
that this classical computation has the interpretation of computing an approximate density 
matrix in some approximate theory. In cases where the boundary is a true asymptotic 
boundary (such as a locally asymptotically AdS boundary) the situation is well understood. 
This corresponds to computing the entropy of a perfectly well defined density matrix in 
the dual field theory . 

Interestingly, there is a simple conjecture for the final answer. The entropy is also 
given by the area of a special codimension two surface in the bulk of the original (n — 1) 
solution. At this surface the circle shrinks smoothly to zero size. The surface obeys a 
minimal area condition. 



S = - nd n [log Z(n) -nlogZ(l)]\ n=1 
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(1.4) 



From now on, log Z(n) denotes the classical gravity action log Z(n) = —SQ rav of the n th 
solution. This formula was first conjectured by Ryu and Takayanagi in the context of 
the computation of entanglement entropy of conformal field theories with gravity duals 
J7j(see H for a review). Proving their formula amounts to proving the above conjecture, 



as we explain below. Notice that (|1.4|) can be viewed as a statement about classical general 
relativity. It is a relation between the actions for classical solutions that are produced by the 
replica trick and the area of the minimal area solution with n = 1. Of course, for solutions 



with a U(l) symmetry, (|1.4|) reduces to the standard Gibbons-Hawking computation. In 
that case, the U(l) symmetry also ensures that the horizon is a minimal surface, with zero 
extrinsic curvature. 

In this paper we will give an argument for ( |1.4| ) based on reasonable assumptions 
regarding the analytic continuation of the solutions away from integer values of n. 



We will also explain why proving ([O]) is equivalent to proving the Ryu Takayanagi 
conjecture. The Ryu- Takayanagi conjecture for the case of asymptotically AdSs pure 
gravity was proven in |)]||. Previous arguments include 0, whose assumptions were 
criticized in [jrU|B 

This paper is organized as follows. In section 2 we perform some explicit computations 
in a simple example. In section 3 we review the derivation of the entropy formula for 
the case with a U(l) symmetry. In section 4 we present the arguments for the main 
formula (|1.4|) . There we explain how the solution looks for n close to one. We also 



derive the minimal area condition for the surface. In section 5 we discuss the connection 
to entanglement entropy in field theories with gravity dual. In section 6 we present the 
conclusions. In the appendices we present some further explicit examples and more details 
on the computations. 

2. A simple example without a U(l) symmetry 

Since our discussion has been a bit abstract, let us discuss a very simple concrete 
example. This example will also motivate some assumptions that we will make later. 



Fursaev Q took the solution for n = 1 and set t ^ r + 2itn everywhere in the bulk. This 
introduces a conical singularity in the bulk. As noted by Headrick flC|| , for integer n, one should 
instead consider solutions which are non-singular in the bulk. 



Let us start with the BTZ geometry 
ds< ^ 



1 + r , 



+ r 2 dr 2 + (1 + r^dx' 



(2.1) 



This metric has a U(l) isometry along the circle labeled by r, r ~ r + 2n. All functions will 
be invariant under translations in x. This direction will not play any role in this discussion 
and we take it to be compact of size L x . Computing the entropy for this solution gives the 
standard area formula, Sq, for this solution. 

We now add a complex, minimally coupled, massless scalar field (j). We set boundary 
conditions that are not U(l) invariant 
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We now compute the gravitational action to second order in r\ for the family of solu- 
tions described above. The metric is changed at order ?y 2 , but since the original background 
obeys Einstein's equations, there is no contribution from the gravitational term to order 
r/ 2 . So, to this order, the whole contribution comes from the scalar field term in the action. 

Namely, for the n th case, we need to consider a spacetime with the same boundary 
conditions as in (2.2) but where t ~ T + 2nn. This implies that the spacetime in the 
interior is 



ds z 



dr' 



+ r 2 dr 2 + (n 2 + r 2 )dx* 



(2.3) 



(n 2 + r 2 ) 

And we need to consider a scalar field in this spacetime. We can write the wave equation. 
The solution of the wave equation that is regular at the origin and obeys ( |2.2[ ) at infinity 

is 



(f> = Ve lT fn{r) 



n n 



fjr) = (nr) n ^ '- 2 F 1 ( -, - H- l;n + 1; 

J v ' K J r (n + 1) V2 2 



nr) 2 ) (2.4) 



Note that f n — ¥ 1 as r — > oo. 

We now evaluate the gravitational action for every n. We evaluate it to second order 
in t/, so we consider the quadratic action for the field <p. Using standard formulas we can 
write 



log Z(n) 



AdS 3 



\V(/)\ 2 = -(2Tvn)L x [r 3 (f>*d 7 



(2-5) 



= (2ttL x ) [1 — nlogn + mfj(n/2) + (linear in n)] 



where L x is the length of the x direction and ijj is the Euler i/j function. The terms linear 
in n include divergent terms that should be subtracted. However, they do not contribute 
to the entropy (|1.3[) . 

We analytically continue in n and compute the entropy via ([T|) to find 



S = S + r,\L x {A-^-) (2.6) 

We can now compare this with the answer we expect from the area formula. This non- 
zero configuration for the scalar field changes the geometry to second order in r\. Thus it 
produces a second order change in the area of the horizon. This change can be computed 
from Einstein's equation. We obtain the same answer (|2.6|) . This is done in detail in 
appendix A, where we also consider a scalar field with an arbitrary mass. 

So, we have explicitly checked the conjecture for this special case. Now, let us make 
some remarks. 

We considered a complex scalar field, but the computation can be done also for a real 
scalar field with boundary conditions <fi = rj cost at infinity. The result is essentially the 
same. See appendix A. 

Notice that the solution for the n th case has a Z n symmetry. This is a replica symmetry 
of the boundary conditions which extends to the bulk solution. So, in this case we are not 
breaking the replica symmetry. Notice that r = is a fixed point of the action of the Z n 
replica symmetry for all n > 1. In this case, the metric has a (7(1) symmetry. However, 
the full scalar field configuration is only symmetric under the Z n H 

Here we have computed log Z(n) and then analytically continued the answer. The 
geometry ( |2.3| ) is well defined also for non-integer n and we can trivially continue it to non- 
integer values of n, and it remains smooth. We could ask whether we can also analytically 
continue the whole field configuration to non integer values of n. Notice that as we vary 
n, the t dependence at the boundary is kept fixed. Thus, even for non-integer n, we will 
keep the same boundary condition. This boundary condition is not compatible with a 
non-integer period for r. We will ignore this. In other words, we will integrate r between 
[0, 2tt] and multiply the result by n. However, as we go to r = 0, we find that the scalar 
field behaves as (p ~ r n e lT , which leads to a singularity for the scalar field at r = 0. The 



In this case there is a £7(1) symmetry which is shift in r combined with a phase rotation 
of the complex field. But for a similar computation with a real scalar field we only have the Z n 
symmetry. 



scalar field, or its stress tensor do not diverge if n > 1. In other words, this appears to be a 
relatively harmless integrable singularity. This singularity seems physically questionable. 
But we are not trying to give a physical interpretation to the solution with non-integer 
n. We are only trying to define it mathematically, as an intermediate step in computing 
the replica trick answer. One could worry that if we allow singularities, then the solution 
will not be uniquely defined. However, we are allowing a very specific behavior which 
determines a unique solution for given boundary conditions. More explicitly, note that 
when we solve the wave equation near r ~ 0, we get two solutions r n e lT and r~ n e lT . We 
set to zero the coefficient of the second solution at the origin. This prescription uniquely 
selects a solution, both for integer and non-integer n. 

This gravity theory in AdS% with a massless scalar field can arise from a Kaluza 
Klein reduction of a higher dimensional theory. For example, it can come from a ten 
dimensional solution of the form AdSs x S 3 x T 4 . Then the massless field can be an 



off-diagonal component of the metric on of the four torus [11,12]. More explicitly, we can 
deform the metric of the four torus as 

ds 2 TA = e 2(t>1 dyl + e 2 ^ 2 dy 2 + e~ 2 ^dy\ + e'^dyj (2.7) 

where <p — <fii + i<p2- We see that the singularity of the field <p at the origin translates into 
a singularity for some of the Riemann tensor components. Let us consider n = 1 + e. Then 
since <p oc r 1+£ e lT this leads to a singularity in some of the Riemann tensor components 
Rai/3i ~ ~ ( no sum over i)i where i are the directions on the four torus and a denotes the 
directions along the two transverse components (labeled also by r and r) Despite these 
singularities the action is finite, as we saw when we computed it explicitly. 

An alternative way to view the solution labeled by n is the following. We consider the 
r circle to have period 2n but introduce a conical singularity at the origin with opening 
angle 2n/n. This is not the same as the gravity solution with n = 1 since the field 
configuration has to adjust to the presence of the conical singularity. Then, when we 
evaluate the gravitational action, we integrate r over [0, 2-7r] but multiply the resulting 
answer by a factor of n. This factor of n arises because the real period of r is 2nn instead 
of 2tc. It is important that we evaluate the gravitational action without introducing any 
contributions from the tip of the conical singularity, since the full space (with the right 
period for r) is non-singular B . This picture makes sense both for integer or non-integer 
n. 



This looks superficially similar to what was discussed in P], but it is different in detail. 
This point of view was also suggested to us by T. Faulkner. 



3. Computation of gravitational entropy when there is a U(l) symmetry 



In this section we will describe the computation of the entropy using Euclidean meth- 
ods in a way that it emphasizes the fact that the contribution comes form the horizon. 
This has been discussed by various authors in a similar form [B JI^ , H[ , ^B JIB , IT] . Here we 
say it two ways that we particularly liked. 



3.1. Entropy from rounded off cones 





+ 
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Fig. 3: A particular combination of geometries that is useful for computing the 
entropy. The first geometry is the correct solution with period 2iin. The last 
geometry contains a conical singularity. It is the solution with n = 1 but with 
the circle identified after t ^- t + 2-rm. For n = 1 the deficit angle of the cone 
is very small and it has been greatly exaggerated here for artistic reasons. The 
two middle ones are identical and correspond to a regularized version of the last 
solution. They only differ for r < a, where a is small regulator. This is not a 
solution, it is an off-shell configuration. All of the configurations obey the same 
boundary conditions at infinity. 



Setting the period of the circle to be r ~ r + 2nn, then we find that the formula for 
the entropy can be written as 



S = -nd n [log Z(n) - nlogZ(l)] n=1 



(3.1) 



Let us consider this expression for n close to one. We interpret the first term in the square 
brackets as the correct, smooth solution when n is not one. We interpret the second term 
as the solution for n = 1 but with a r which has period t ~ t + 2nn. This solution has a 
conical singularity at the origin. However, we do not include any contribution from the 



conical singularity. We simply integrate the gravitational action density away from the 
tip. 

We now evaluate the difference in the square brackets in ( |3.1|) by adding and subtract- 
ing a smooth geometry which is the same as that of the cone far away from the origin, but 
it is a regularized cone near the origin, see fig. 3. This smooth geometry is not a solution 
of the equations of motion of the theory, it is an off-shell configuration. We are simply 
introducing it to help us perform the computation. It is possible to choose this off-shell 
configuration in such a way that the metric differs only by an amount of order n — 1 from 
the true solution. 

Thus we get 

S = -nd n [(log Z(n) - log Z oS (n)) + (log Z oS (n) - nlogZ(l))] n=1 (3.2) 

Each of the terms in the brackets is the action for one of the configurations in fig. 3. Since 
the off shell configuration that corresponds to a regularized cone differs by a first order 
term in n — 1 from a solution of the equations, we see that we can interpret the first 
parenthesis as the result of doing a first order variation away from a solution (the solution 
with period n). This first order variation vanishes due to the equations of motion for the 
solution with period n. Notice that both metrics obey the same boundary conditions at 
the boundary, so that there are no boundary termsH 

So all that remains is the second parenthesis. The second parenthesis contains the 
difference between a smooth cone and a regularized cone. This receives a contribution only 
from the region near the tip of the cone. This contribution is extensive in the area of the 
horizon, namely the area of the surface transverse to the tip of the cone. The region near 
the rounded tip of the cone contains an integral of J d 2 x^fgR along the cone directions 
which gives 

d 2 x^gR ~ 4tt(1 - n) (3.3) 



JRi 



Reg Cone 

Thus, the final answer has the form 



■• , . - ( ~ f ,2 — ~"\ Area 



S = —-—(Area) -nd n / d'xy/gR = — (3.4) 

i07rCrAr \ J Reg Cone J 4Grjv 



The absence of boundary terms is clearest if we write the action in a non-manifestly covariant 
form using only first derivatives of the metric. Then the fact that the two configurations obey the 
same boundary conditions for the metric implies that there are no boundary terms. 
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One can consider a metric that explicitly regularizes the cone, such as [|16 



ds 2 = dr 2 g 2 (r)+r 2 dr 2 (3.5) 

where g = n + o(r 2 ) at r ~ and g = 1 for r > a, where a is a small distance which sets 
the size of the regularization. Inserting this metric into the gravitational action we get 
(|3.3|) . One can choose a completely explicit function such as g = 1 + (n — l)e~ r l a , for 
example. In this case we can see explicitly that the metric perturbation is of order (n — 1). 

3.2. Entropy from apparent conical singularities 

Another way to think about this problem is as follows. First we note that, since the 
solutions are invariant under time translation, the evaluation of log Z{n) is the same as 

logZ(n) = n[logZ(n)} 2n (3.6) 

where [logZ^n)^ is the gravitational action density for the solution labeled by n but 
integrated over r from [0, 2tc] (instead of [0, 27m]). We can now write the entropy as 

S = -n 2 d n [logZ(n)] 27T (3.7) 

Note that the solution labeled by n is a smooth geometry if the r circle has period lixn. 
On the other hand, imagine we wanted to view it as a configuration where the r period 
continues to be 2-k. In that case, it is a geometry with a conical singularity whose opening 
angle is 2n/n. Thus we can view 

[log Z(n)] 27r (3.8) 

as the gravitational action of a configuration with r = r + 2-k but with a conical singularity 
with opening angle 27r/n, without including any curvature contribution from the conical 



singularity. Then we see that the expression of the entropy (|3.7|) involves taking a derivative 
with respect to n. When we change n we are changing the opening angle of the singularity. 
In addition, we are changing the metric and other fields everywhere since they have to 
adjust to this new strength of the singularity. However, since the original solution (the 
solution with n = 1) is a solution of the equations, we would naively expect that a first 
order variation of the metric and other fields should vanish due to their equations of motion. 
This naive expectation is essentially right, except for the fact that we are changing the 
boundary conditions at the origin, since the strength of the conical singularity is being 
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changed. Thus, the only change in the action comes from a boundary term. In other 
words, when we change n the action changes as 



-d n [\ogZ(n)} 2 ir\n=i = / E 9 d n g + E<j>dn(f>+ 

+ ^~ I dy D - 2 ^(Wd n g^ r - g^Vrdng^ 



i r a (3 - 9) 
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where E g and E§ are the equations of motion for the metric and other fields, which vanish. 
Here y are the coordinates along the r = surface. The boundary term vanishes at the large 
r boundary since we are choosing boundary conditions in such a way that the variation 
of the action gives the equations of motion without extra boundary terms. On the other 
hand, at the horizon (at r = 0), we do get a contribution from the boundary term. This 
boundary term produces the area contribution. Note that the n derivatives of the metric 
are evaluated at the horizon. For example, in the parametrization ds 2 = n 2 dr 2 + r 2 dr 2 
near the origin, we get, as the only non- vanishing component, d n g rr \ n= \ = 2. With these 
expressions we can evaluate the parenthesis in ( |3.9| ) and obtain 2/r. 

This derivation easily generalizes to theories with higher derivative actions, giving the 
Wald entropy H^g. 



Note that in both cases we used explicitly the locality of the action along the r 
direction. It would be interesting to find the corresponding formula in weakly coupled 
string theory exactly in a' . 



4. Argument for the entropy formula (1.4) 
4-1. Properties of the metric for n integer 

For n = 1, the boundary contains a circle which we label by the coordinate r. Recall 
that the boundary is the surface where we are putting boundary conditions. This circle is 
non-contractible on the boundary, but it can be contractible in the interior of the geometry. 
Here by boundary, we mean the boundary where we set boundary conditions for the 
gravitational action. It need not be an asymptotic boundary. 

The metric and all fields are periodic on this circle. Let us collectively denote these 
fields as i/j(t), with 

i/j(t) ~ ip(r + 2tt) (4.1) 
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Of course, the fields depend on other coordinates, but here we are highlighting their r 
dependence. We impose boundary conditions 

^(t) boundary = ^b{t) , V>b(t) = $b{t + 2%) (4.2) 

where we specify the functions i/jb(t), which are periodic. 

The solution with n > 1, has exactly the same boundary conditions (|4.2|) , but we 
require the periodicity t = r + 2-kti on the r circle. This implies that the boundary 
conditions have a Z n symmetry. We assume that the bulk solution continues to have this 
Z n symmetry.^. 

Each of the solutions for n > 1 has a special co dimension two surface which is left 
invariant by the action of Z n . We will focus on this surface. We can choose a coordinate r 
which is a radial coordinate away from this surface and an angle r. The true angle around 
the surface is really a = r/n, we have chosen r to have the same period as the one we 
have at infinity (a ~ a + 2ir). The metric in the two directions transverse to this surface 
has the form 

ds 2 = n 2 dr 2 + r 2 dr 2 + ■■■ (4.3) 

where the factor of n comes from demanding that there is no singularity at r = 0. In 
addition, all fields are required to have an e lkT dependence, with integer k. This comes from 
the period of r and the Z n symmetry. Thus, a scalar field would behave as r n e lT ~ r n e ina 
near the origin, as results from demanding that it is non-singular. 



There are cases where the Z n symmetry can be broken. This was observed computing 



the Renyi entropies of circular regions in the three dimensional interacting O(N) model [21]. 
We thank I. Klebanov for pointing this out. The simplest gravity solutions can also develop 
instabilities. For example, if one considers gravity in AdSd+i with a boundary Hd-i xS 1 . If 
the radius of S 1 is equal to the radius of Hd-\ then the full solution is AdS, viewed as a black 
brane with a hyperbolic spatial section. If we make the S 1 n times larger, then for large n, we 
approach an extremal black hole with an AdS2 X Hd—i near horizon geometry. This can lead to 
bad tachyons for rn 2 R\ dS < —d/4. Thus if the original AdSd+i has tachyons in the allowed 
range — d 2 /4 < m 2 E? AdS < — d/4, then we will have an instability. This is similar to the discussion 



of [22], where an extremal Reissner-Nordstrom black brane in AdS^, with a near horizon geometry 
AdS2 x R 2 was considered. Good tachyons in AdS^ can be bad tachyons in the AdS2 region if 
— | < m 2 R 2 AdS < — | . However, in suitably simple cases we expect this replica symmetry to be 
present, and our arguments will apply. 
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As a side remark, notice that if the bulk space has no fixed points under the Z n action, 
then this means that we can choose the coordinate r in the interior so that this circle never 
shrinks. An example is a space with topology K d ~ x x S 1 , but with a metric that depends 
on the coordinate along the S . In these cases the entropy is zero. The reason is very 
simple, the solution for the n th replica is the same as the solution with n = 1 but with 
a longer circle so that log Z(n) = nlogZ(l). Here, of course, we used the locality of the 
classical action. 

4-2. Metric for n non-integer 

Here we make some assumptions on the form of the metric when n is not an integer. 
We will continue to impose exactly the same boundary condition Q4.2|) , which is periodic 
with period 2tt. This is not compatible with t ^- t + 2nn. Now, in the region where the r 
circle has positive size, we can ignore this problem and think of r as being non-compact. 
When we evaluate the action, we can integrate the r direction from to 2n and then 
multiply by n. 

However, we expect that there is still a surface where the r circle shrinks to zero. 
For the two dimensions transverse to this surface we impose that the metric continues to 
behave as in Q4.3| ), even though n is not an integer. The rest of the fields, including other 
components of the metric, are chosen so that they are periodic in r — > r + 2n as in ( f4.1| ). 
This implies that the field configuration is singular at r = 0. However, we expect that this 
singularity is as harmless as the one we had for the scalar field in section 2 . 

This is seems a reasonable assumption. As evidence for its validity we can point to 
the explicit example mentioned in section 2 . 

An equivalent way to specify the solutions is to compactify the r circle to r + 27r in 
all cases (all values of n) and demand that there is a conical defect angle with opening 
angle 2n/n in the interior. We do not introduce any contribution to the action from the 
tip of the cone. In addition, we multiply the gravitational action by a factor of n. This 
is mathematically equivalent to what was discussed above and the reader can choose the 
preferred interpretation. 

Note that this is similar to introducing a cosmic string (or cosmic D — 3 brane) with 
opening angle 27r/n in the original solution, with the metric backreacting as necessary to 
account for its presence. 

As n — > 1 the solution goes over to the solution with n = 1. Thus, this analytically 
continued solution is close to the n = 1 solution and we can expand it in powers of n — 1. 
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4-3. Derivation of the minimal surface condition 

We emphasized that for n > 1 we have a special surface where the circle shrinks, 
and is fixed under the Z n action. But for n = 1 there is no obvious special surface, 
since there is no unique way to choose the coordinate r in the interior once it is not 
associated to a U(l) isometry. So, when we expand the solution in n — 1 we need to select 
a surface. Motivated by the Ryu-Takayanagi conjecture we want to select a minimal area 
surface. In this subsection we will explain the origin of this minimal area condition. The 
final conclusion is that the condition comes from demanding that the solution obeys the 
Einstein equations to leading order in n — 1. This derivation is essentially the same as 
the derivation of the equations of motion for a cosmic string (or D — 3 brane) from the 
behavior of the metric near the conical singularity. This problem was analyzed previously 
^Ell- 



in 



Two dimensional dilation gravity 

It is good to start with a simple situation first. For that purpose we will consider a 
two dimensional dilaton gravity where the action is 

-Scrav = ^ / d 2 x^e~ 2 ^ [R + 4(Vy?) 2 + • • •] (4.4) 

where the dots indicate other fields, or a potential for f, etc. Notice that if we have a 

solution with a horizon, then e~ 2ip at the horizon plays the role of the area in Planck units 

of the higher dimensional gravity solutions. In this case the codimension two surface is just 

a point. The minimum area condition is that e~ 2ip is a minimum (or really an extremum) 

at this point. We will derive this condition from demanding that the configuration for 

small e = n — 1 obeys the linearized field equations near r = 0. In other words, expanding 

the fields around the n = 1 solution, and assuming the periodicity condition for the fields, 

( fl.l|) , we will see that we can only obey the equations if di<p = 0. 

Let us say that as n — > 1 the special surface goes over to some point of the n = 1 

manifold. Let us pick this point to be the origin in some coordinate system x l ,x 2 . Then 

the metric of the n = 1 solution around this point is ds 2 = dx\-\-dx 2 ,-\-o{x 2 ). The field (p is 

regular at this point. Now, for n— 1 = e we expect a metric of the form ds 2 = e 2p {dx\-\-dx\), 

with e 2p = r^™ -1 ), as r — > 0. Then to first order in e we have 5p ~ — elogr to be the 

first order solution. We consider the two following equations for two dimensional dilaton 

gravity 

=e~ 2if (4d^d z p + 2d%<p) + T z 



-imatter 
■ zz 



= e - 2tp {4d^d sP + 2d\y) + T, 
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matter 

zz 



(4.5) 



where z = x 1 + ix 2 , z = x 1 — ix 2 . Here 7" matter denotes the stress tensor for the rest of the 
fields of the theory, coming from the dots in (fOJ) . Expanding the first equation to first 
order we find 

-2d,v?(0) - + 1d\b<p + 5T™ atter = (4.6) 

and a similar equation by expanding the second. Here d z <p(Q) is the derivative of the field 
for the n = 1 solution at the origin. It is just a z independent constant. Since the matter 
stress tensor is not expected to be singular at order 1/r, we find that 

2 ed z <p(0) 2 ed s <p(0) 
a z cnp oc , o 2 0ip oc = (4.7) 

z z 

up to terms that are less singular as r — > 0. Now we assume that the solution for Sip has 
a fourier expansion with integer powers of e lT . The first equation in ( |4.7| ) suggests that 
we try a solution proportional to Sip oc zlogz. However, the periodicity condition under 
shifts of r suggests that we should consider Sip oc zlog(zz). However, this is not a solution 
of the second equation. Thus this implies that the gradients of the field should vanish at 
the origin. 

More formally, we can argue as follows. The periodicity condition implies that if we 
take the r derivative of any field and integrate over r between zero and 27r, we should get 
zero. This is true both for 5<p and its derivatives. In particular, note that the following 
combination of derivatives gives 

d T [(rd r — l)d z 8ip) oc (zd z — zdz)(zd z + zd z — l)d z 5ip oc ed z ip(0) (4.8) 



where we used both equations in ( f4.7|) . Now the integral over r of (|4.8| ) should be zero, 
according to our assumption about the periodicity of Sep. This then implies that d z ip(0) = 
0. 

In summary, in this case we found that the condition comes from the zz and zz 
components of the Einstein equations. In higher dimensions we expect that this will come 
from Einstein's equations in the directions normal to the surface. 

Note that if we changed the coefficient of the dilaton kinetic term in (|4.4|) from 4(V<p) 2 
to (4 + cr)(Vc£>) 2 , then we would be adding terms of the form ad z ipd z ip to the equations in 
( |4.5|) . Expanding around the background solution such terms lead to contributions that 
are subleading, in the expansion around the origin, compared to the terms already taken 
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into account in (4"7F). Thus, if we had a two dimensional action with a different coefficient 



for the dilaton kinetic term, we would have reached the same conclusions. 

Einstein gravity in D dimensions 

We now go back to the case of Einstein gravity. In general, we can expand the metric 
of the n = 1 solution around the special surface as 

ds 2 =dr 2 + R 2 dr 2 + bidrdy 1 + g i jdy' l dy j , 

9ij = hij + r cos tK}j + r sin rKfj + o(r 2 ) (4.9) 

R=r + o(r 3 ) , bi = o(r 2 ) 

where r is coordinate normal to the surface and yi are coordinates along the surface. Here 
Kfj are the two extrinsic curvature tensors, hij depends only on yi but not on r or r. 
When we deform away from n = 1 we assume that we cannot change the period of the 
cosines above. 

When n = 1 + e some of the metric components generically go like r 1+e . This can 
give rise to terms in the equations of motion going like 1/r. These terms can only come 
from situations where we have two derivatives along the transverse directions (the r and 
t directions) . Such terms in the equations of motion are the same as the ones we would 
obtain by performing a dimensional reduction from D dimensions to the two transverse 
directions. This brings us back to the previous case. More explicitly, we write the full D 
dimensional metric as 

ds 2 = e 2p (dxj + dx 2 2 ) + e~^g ij dy i dy j + o(r 2 ) , det(^) = 1 (4.10) 



where c/ij is the transverse metric appearing in ( f4.9| ) but rescaled so that its determinant 
is one. The off diagonal terms in ( ^.9| ) do not contribute to terms of order 1/r in the 
equations of motion. Both g^ and (p depend on all the coordinates, the yi and the Xi. Here 
we have pulled out the overall volume factor of the transverse space and parametrized it by 
f. Dimensionally reducing to the first two dimensions gives us ( fl.4| ), but with a different 
coefficient for the dilaton kinetic term. Thus, we obtain the same conditions that <9 X « <p = 



for a = 1,2. Now if we translate between cp and the original metric (|4.9|) we find that 



-4<p = log(det(hij)) + x x K x + x 2 K 2 + o{r 2 ) , K a = h ij K% (4.11) 



This is to be expected since this coefficient can be changed by a field redefinition of the 
metric. 
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where K a are the traces of the extrinsic curvature tensors. We then see that the condition 
d a ip = implies that 

K 1 = K 2 = (4.12) 

Namely, the traces of the extrinsic curvatures should vanish. There are two directions 
that are transverse to the surface so we have two relevant extrinsic curvatures. These 
coincide with the equations of motion for a minimal area surface. There are two transverse 
directions to the surface and thus two equations. In appendix B we derive ( |112|) directly 
in D dimensions, without doing the dimensional reduction. 

Note that the non-trace part of the extrinsic curvatures are not constrained to vanish. 
In fact, already in our simple example of section 2 we have non- vanishing extrinsic curvature 
if we interpret the scalar field as coming from a component of a higher dimensional metric 
as in O-0 

4-4- Computation of the entropy using the cone method 

Once we have established the form of the solution, we can compute the entropy using 
the cone method as explained in section 3. The arguments are similar, but one has to 
check that the mild singularities we discussed above cause no problems. 

Let us discuss this first for the case of AdS plus a scalar field discussed in section 
2. There, the singularity is only present in the scalar field which behaves as </> ~ r c re lT . 
With this mild singularity, if we integrate by parts in order to use the equation of motion 
for </>, it is clear that we do not run into any problem at r = 0. The most dangerous 
term seems to come from the variation of 5 f drr(d r <fi) 2 = 2 J drrd r (j)d r b~(f) — > rd r <pd(f)\o- 
However, 5<p would also vanish at the origin if we are considering the variations that come 
from varying n in the solution. In other words, when we compare the correct configuration 
with n — 1 > and a regularized cone, we can consider a regularized cone where <p has 
the same type of singularity at the origin. This shows that the first parenthesis in fl3.2| ) 
vanishes. 

The second parenthesis only gives us something interesting if we consider terms that 
have two derivatives acting on the metric, otherwise their contribution is going to be small 
as we remove the regulator. Thus, only the metric in the two directions transverse to the 



More explicitly, in the notation of that section, if the field at the origin goes as (p = re %T this 
leads to cf>\ = rcosr = x 1 and there is an extrinsic curvature component K 1 1 ± = —K 1 3 3 using 



the coordinates fl2.7| ) and (|2.1| ) for the 3-d part. 
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minimal surface are relevant. And in those dimensions the computation reduces to the 
usual one, with the contribution coming only from the curvature term in the action. 

In conclusion, evaluating the differences in ( |3.2| ) we find that the answer is equal to 
the area, as we wanted to prove to argue (|1.4|). 

The discussion so far was completely local in the directions transverse to the "horizon" . 
Here by horizon we mean the point in the two transverse directions where the circle is 
shrinking to zero size. In some cases this "horizon" can have multiple disconnected regions. 
Then, we should sum over the areas of each of the horizons. Even when we have multiple 
horizons, the period of the r circle is the same in the whole solution. 

4-5. A comments on other U(l) symmetries 

Throughout this discussion we have focused on the particular geometric circle that 
we used to define the density matrix and the replica trick. We considered cases where 
we have no translation symmetry along the circle. However, we can have other U(l) 
symmetries. As a simple example, we can have a U(l) gauge field in the bulk. Then as in 
the ordinary case, 0] , the integral of the gauge field along the r circle should vanish at the 

origin j Q drA T = (up to global gauge transformations) , where the r circle shrinks. 

i — >o 

This should hold for all n, both integer and non-integer. At the boundary we can fix the 
holonomy of A along the r circle as we please. In order to compute the entropy of the 
density matrix with a chemical potential we should fix the integral p = L A T at the large 
r boundary. If we keep everything else fixed at the boundary but we vary p, this has the 
interpretation of changing the density matrix p —¥ e 1 ^ p where Q is the charge associated 
to the U(l) symmetry. We can compute the entropy of this density matrix by treating 
this boundary condition as we treated all other boundary conditions. Namely, A T (r) is 
kept fixed. Therefore its circulation over the r circle of length 27m is np. Of course, if the 
holonomy in the r circle is different at the origin (r = 0) than at the boundary, then we 
will have a non-zero field strength in the bulk. The computation of the entropy is identical 
to what we discussed in general. 

This other U(l) symmetry can also be an ordinary geometric isometry, and its treat- 
ment is similar. 
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5. Connection to the Ryu-Takayanagi formula 

We presented the computation of the entropy of the gravitational density matrix in a 
form that is very general. The objective was to emphasize that (|1.4| ) is really a statement 
about an analytic continuation of classical solutions. In this section we explain why the 
conjecture ( |1.4| ) for the entropy is related to the Ryu-Takayanagi formula for entanglement 
entropy. 

The Ryu-Takayanagi formula is a conjecture in the AdS/CFT context 0. In the 
quantum field theory one is interested in computing the entanglement entropy of a spatial 
region A on the boundary of the field theory. This spatial region has a boundary d A. The 
conjecture is that this entanglement entropy is given by a the area (in Planck units) of a 
codimension two minimal surface in the bulk whose boundary ends on dA. 




Fig. 4: The Ryu-Takayanagi conjecture. The entanglement entropy of a region 
A in a conformal field theory is given by the area of a minimal surface in the bulk 
of AdS that ends on dA (the boundary of region A) at the boundary of AdS. 



In principle, we can compute the entanglement entropy of the region A by using the 
replica trick [p5| , |26|| . This is a general method for computing entanglement entropy in 
quantum field theories. The idea is to take n copies of the field theory and match them 
together so that by moving in a circle around d A we go from one copy to the next. Going 
n times around this circle we come back to the original copy. Thus at dA there is a conical 
defect with a 2nn opening angle. This appears to be a singular metric. However, one can 
choose a conformal factor that diverges at dA in such a way that the size of the circle 
around dA is finite. 
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a 



b) 



Fig. 5: Geometries that we are considering to compute the entanglement entropy 
in the field theory, a) Semiplane. Region A is half of the plane and its boundary, 
dA, is at xo = xi = 0. b) xo, x\ view of the semiplane and the coordinate r. 



X 



Xq 











/ 


T 




T 


/xi 


C T>^' / 


/X>2 




/ /X\ 




//X\ 





a) b) 

Fig. 6: Other geometries in three dimensions, a) Disk configuration, b) Slightly 
deformed disk. The r coordinate goes around dA. 



This is most easily understood for simple regions [p7| , |28[| . Imagine we have a conformal 
field theory R 1 ' d ~ 1 at the boundary. Then we can choose a region A defined by x\ > 0, see 
figure fig. 5. The boundary of the region is the surface xq = 0, x\ — 0. Going to Euclidean 
space, R d , we can combine the directions x$ and X\ into two directions labeled in polar 
coordinates by r and r. The metric is 



ds' 



r dr + dr + dx 



Ji 



-> dr 1 + 



2 . dr 2 + dx 2 



(5.1) 



where x are the rest of the spatial coordinates. In the right hand side of (|5.1|) we have 
multiplied by an overall conformal factor 1/r 2 to put the metric in the form of S l x H^-i- 
We can now easily perform the replica trick, it corresponds to changing the length of 
5" 1 from 2n to 2nn. Clearly this metric, S^ x H^-i is a perfectly legal metric and we 
can consider its gravity dual. It is a certain black brane. In this case, we have a U(l) 
isometry in the rescaled coordinates and then the entropy computed using the replica trick 
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or using the ordinary Gibbons-Hawking formula is exactly the same. Note that at the AdS 
boundary the circle S 1 has a nonzero size everywhere. In the interior of AdS it shrinks to 
zero at a "horizon" . Notice, in particular, that the half space region we discussed above 
can be conformally mapped to a spherical region Yl x l — 1- I* 1 this case, the circle S 1 
appearing in (|5 . 1 [) corresponds to a coordinates that goes around dA as in figure fig. 6. 



Now, this was a very simple region. If we consider more complicated regions, then 
it is not possible to choose a system of coordinates and a conformal rescaling such that 
the metric is independent of the angular direction r. In all cases we will have an angular 
direction, r, since it is the direction we used to perform the replica trick construction. The 
choice of this coordinate is completely arbitrary, as long as it goes around the boundary 
of region A. As we go near dA we have a problem which locally looks like ( |5.1| ), and we 
can choose a conformal factor which makes the metric non-singular as in ( |5.1| ) for all the 
replicas. The difference with (|5.1| ) is that, as we increase r, we will have extra terms in the 
metric that can have some r dependence. This dependence always involves powers of e ±ir 
since this is just the statement that the r direction is parametrizing circles in the original 
boundary geometry. All the statements in this paragraph involve the boundary geometry, 
the geometry where the field theory lives. These replica trick boundary geometries simply 
amount to letting the circle have size t ~ t + 2nn, without changing any of the functions 
that appear in the boundary geometry. All such functions are periodic under r — > r + 2-7T 
0. 

Thus, the field theory replica trick, can be translated, via the standard AdS/CFT dic- 



tionary |2t|[>0] , to a problem in gravity which is identical to the problem that we discussed 
in section 4. Here no conjecture is involved other than the original AdS/CFT relation. 
The replica trick then defines the entropy as in ( |l.ij| ). In order to do that, we need to 
analytically continue in n to n ~ 1. 

The Ryu-Takayanagi conjecture boils down to a statement in classical geometry. It is 
the statement we discussed in section 4. Computing log Z(n), using smooth geometries, 
analytically continuing in n, and computing the entropy defined in (|1.3| ) gives the area 
formula in ( |1.1| ). 

Notice that in a setup where A is a spatial region contained at x° = on the boundary, 
then there is a time reflection symmetry a; — > —x°, which translates into r — > —r for the 



If we want to explicitly parameterize the metric in this way, we might need to choose different 
coordinate patches, as usual. 
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circle in the Euclidean solution. This implies that we can go to Lorentzian signature, as 
usual with x° — > ix°. This translates into r — > it . Now the region where the r circle 
is shrinking to zero corresponds to a horizon in the bulk. It is a horizon for an observer 
sitting at fixed small r. 

There is a generalization of the Ryu-Takayanagi conjecture for situations that are 
time dependent |H]. It again involves an extremal surface ending on cL4, but in the full 
Lorentzian spacetime. In those cases there is no obvious Euclidean continuation to perform 
the replica trick. This suggests that there should be a way to think about the problem 
which does not go through the Euclidean solutions and the replica trick. We should remark 
that in some cases we can perform a replica trick in the Euclidean geometry for regions that 
depend on the Euclidean time and then one can analytically continue to the Lorentzian 
signature solutions. Some examples were discussed in |32] . 



5.1. General entanglement interpretation 

In the introduction, we presented the computation of the generalized gravitational 
entropy as a property of the density matrix constructed by integrating over a circle in 
Euclidean time. It is natural to ask whether there is a general Lorentzian interpretation 
that involves entanglement. This is indeed the case in the Ryu-Takayanagi discussion of 
entanglement of a subregion of the boundary. 





B 



(a) 

Fig. 7: We consider periodic boundary conditions with a reflection symmetry 
t — > —t. In (a) we see that by cutting at r = we get a density matrix p ac , where 
a and c label the states on the two sides of the cut. In (b) we note that we can 
cut along the moment of time reflection symmetry r = 0, tt. Then we get a pure 
state in two separate Hilbert spaces labeled by A and B. The bottom half of the 
picture can be viewed as a state x/j a b and the top part as ^ cd . Tracing out over the 
B Hilbert space, we recover p ac = Y2h' l ^ ab '^'cb- At this moment of time reflection 
symmetry we can also continue to Lorentzian signature. 
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Here we would like to point out that in very general situations we can also have an 
entanglement interpretation. Suppose that the boundary conditions have a moment of 
time reflection symmetry Say that this acts as r — > —r. Then by cutting the boundary 
conditions at r = 0, 7r we can interpret the lower part of the evolution as specifying a pure 
state |\P) in the product of two theories, which we call A and B. See fig. 7. Similarly, 
the upper part can be viewed as specifying the state (\I/|. The density matrix can then 
arise by tracing over one or the other subsystem. And the entropy can be interpreted as 
entanglement entropy for system A with B. This is the same as in the eternal black hole 
discussion |33|,[34 . 




horizon 



r=0 

A ■ 

Fig. 8: Here we consider a situation with asymptotically AdS boundary condi- 
tions. The boundary conditions contain a small time dependent deformation which 
vanishes at infinity. So in the far future we settle down into a stationary black hole 
on both sides. The entropy of these black holes is bigger than the entropy of the 
initial entanglement since, the time dependent boundary conditions have sent in 
energy and have increased the entropy of the system. In other words, there was a 
non-zero flux of energy through the horizon which increased its area. The dotted 
lines indicate the matter falling through the horizon. 



The bulk solution is also expected to have a time reflection symmetry in this case. 
Under r — > it we get a Lorentzian solution. The vicinity of r = looks locally like 
Rindler space. This procedure generically produces a time dependent solution and we 
might get singularities or horizons in the boundary conditions. We can consider a situation 
where the Lorentzian time evolution can be performed out to infinite time without ever 
connecting again the two boundary regions or encountering singularities on the boundary. 
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An example is the following. We start from a Euclidean black hole but with a small 
perturbation of the boundary conditions which is smooth in Euclidean time and goes to 
zero at large Lorentzian time. More concretely, we can consider the model of section 2 and 
set the boundary conditions 6b = ox COST ' , When we go to Lorentzian time this becomes 
^ B = 2+cosh 2t wriere the ± corresponds to the A and B sides respectively. Note that 
these go to zero at large times. We expect that solution should be qualitatively like fig. 8. 
In cases that arise from entanglement of subregions via AdS/CFT, the fact that the 
causal horizon is closer to the boundary than the minimal surface that computes the 
entanglement entropy was noted in |JT|(see also |]35| , p6| , |37| j3! 



6. Conclusions and discussion 

In this article we have noted that we can generalize the concept of Euclidean gravita- 
tional entropy to more general situations than the ones associated to thermal equilibrium. 
In particular, we have considered euclidean solutions that contain a circle t ^- t + 2n. 
We have introduced a boundary, setting boundary conditions which are r dependent but 
periodic under t ^- t + 2n. Thinking of gravity as a holographic theory, we view these 
boundary conditions as defining the system. Euclidean evolution on the circle produces 
an un-normalized density matrix. The Euclidean solution gives us the trace of this density 
matrix. By performing the gravity version of the replica trick we have defined traces of n th 
powers of the density matrix. These are geometries with exactly the same boundary con- 
dition as functions of r, but where the r variable is taken to have period t ^- t + 2nn. For 
integer n the bulk geometries are smooth and free of any conical defects. These geometries 
are computing the trace of the n th power of the density matrix. By analytically continuing 
in n and taking a derivative near n = 1 we can compute a quantity that is interpreted as 
the Von-Neumann entropy of the underlying density matrix. Note that all computations 
are classical. The density matrix we are talking about is a hypothetical density matrix 
in some underlying theory of quantum gravity. In AdS/CFT situations we can give an 
precise definition for this density matrix. 

A version of the Ryu-Takayanagi conjecture is that this generalized gravitational en- 
tropy, computed in this fashion, is given by the area of a minimal area surface in the 
original geometry (the solution with n — 1). 

We have given some arguments for the correctness of the Ryu-Takayanagi conjecture. 
The arguments involved the assumption that we can analytically continue the geometries 
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away from integer values of n. We further made the assumption that these analytically 
continued geometries, for small e = n — 1, are smooth in the two directions transverse to the 
minimal area surface but can have mild singularities which are not important for evaluating 
the action. We do not view these metrics as physically meaningful, we view them just as 
a tool for deriving the Ryu-Takayanagi formula. Our assumptions were motivated by 
considering a simple example, described in section 2. But we have no further justification 
other than the fact that they hold in this example and seem reasonable assumptions. We 
have derived the minimal area condition by demanding the existence of a small deviation 
away from the n = 1 solution that is consistent with our assumptions on the type of 
singularities that are allowed. One simple way to state the type of allowed singularities is 
to do a dimensional reduction of the whole configuration to the two dimensions transverse 
to the minimal surface. Then we have a two dimensional metric, a dilaton field that 
multiplies the two dimensional curvature in the action and a set of other fields. Then the 
metric should be smooth and the gradient of the dilaton at the minimal surface should be 
zero, which is the minimal area condition. All other fields can have mild singularities of 
the form (j) ~ z\z\ at the origin. When n is not an integer we evaluate the gravitational 
action by integrating r between [0, 2n] and then multiplying by n. We have also argued 
that this method gives rise to the area formula for the entropy, essentially for the same 
reasons as for the case with the U(l) symmetry. One way to understand this is that all 
non £7(1) invariant fields are going to zero at the origin. Then the methods described in 
section 3 give the usual formula. 

An alternative way to view the solutions is to imagine that we keep the original period 
of the circle, r ~ r + 2n but we introduce a cosmic string (or cosmic D — 3 brane) with 
a 27r/n opening angle. In addition, we multiply the resulting action by a factor of n. For 
n close to one we have a very light cosmic string that deforms the geometry very slightly. 
We can then view the entropy formula as arising from the Nambu action for this cosmic 
string. Also the minimal area condition comes from minimizing this Nambu action. The 
long and detailed discussion that we presented tried to justify these statements in detail. 

One interesting open question is whether one can generalize the derivation to the 
time dependent case considered in [[JTJ, where, generically, there is no obvious Euclidean 
continuation. 

Another interesting direction is to generalize the discussion to gravity with higher 
derivatives. The most naive conjecture is that the entropy is given by the Wald formula. 



However, this conjecture was argued to be wrong in |S9[, where a modified conjecture was 
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made for the case of Lovelock gravity. A more informed conjecture is to say that we get 



the Wald-Iyer formula proposed in section 7 of |]l9|j . In fact, this reduces to the proposal 

in [[J9| for Lovelock gravity. It would be interesting to see whether this is correct and what 

the equations for the surface are. 
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Appendix A. Example of a scalar field in AdSs 

In this appendix, we consider a massive scalar field in AdSs and show explicitly that 
the entropy that we compute using the replica trick is equal to the modification of the area 
due to the presence of a non-zero scalar field background. 

A.l. Massive scalar field 

For a massive scalar field we have equations which are very similar to the ones in the 
text. We consider a complex scalar field of mass m. Setting the radius of AdSs to one we 
need to impose the boundary condition 

0| rc = V e iT r^- 2 , A = 1 + y/m? + l (A.l) 

where A is the scaling dimension of the corresponding operator and r c is a large value of 
r which represents the cutoff surface. The relevant solution of the wave equation on the 
metric ( |2.3|) is 

4, = ^^-rt 2 , /(r) = r" 2F1 (l ~ £ + 1, 1 + £; n + 1; -r 2 ) (A.2) 



f(nr c ) c '^'-^22 '22 

We then evaluate 

logZ(n) = - Jd 3 x^[\VcP\ 2 + m 2 \ ( f>\ 2 ] = -(2ivn)L x ^ r r\d r 
=(2nL x )\r]\ 2 [B(n,A) -flinear in n] 



(A.3) 



where the terms linear in n also include all divergent terms. It is important that these 
counterterms do not give rise to any non-trivial n dependence. This is due to the fact that 
we keep the r-dependence of the boundary conditions fixed as we vary n. We also defined 
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, AN 2n 3 " 2A r 2 - A r (n±A 
S(n, A = A V 2 ; A.4 



We can then compute the entropy to order n 2 from ( |1.3|) , which gives 

S\ v 2 =-nd n [logZ(n) -nlogZ(l)]\ n=1 = 

2 f47r[2(A-2)A + (l-A)7rtan(7rA/2)]r(2-A)r(^±i) 2 l (A.5) 

~~ V I r(^) 2 r(A) J 

AJ2. Change in the metric from Einstein's equations 

Now we will study the backreaction of the scalar in the metric. The action is 

S= [ [l?-2A-|V</>| 2 + m 2 H 2 ] (A.6) 

JAdS 3 

with A = — 1. The equations of motion are 

R^- 9 -f{R-2)=T {ilv) (A.7) 

where T^„ = d^cffd^cf) — ^p(|V</>| 2 + m 2 \4>\ 2 ). The ansatz for the metric is 

ds 2 = —,, —, -dr 2 + (r 2 + l) (1 + v(r))dx 2 + r 2 dt 2 (A.8) 

where g(r),v(r) will be 0{n 2 ). If we expand Einstein equations to first order we obtain 

three equations for the diagonal components. There are only two independent equations 

since the last one will give us the scalar wave equation when the first two are satisfied 

2r 



g'{r) = T x 



(r 2 + 1) 

2rg(r) ' 



v'(r) = 2rT rr — ., 

V ' (r 2 + l) 2 

Since we consider a configuration with d x (p = 0, we can relate the components of the 

stress energy tensor: T rr = (d r (j)) 2 + (1 + r 2 )~ 2 T xx . We then find 

v'{r) = 2r(d r <P) 2 + |- (^;) (A10) 

And 

/"OO /"OO 

u(0) = -2/ drr\d r (p\ 2 -+ S\ v 2 = 4n5A = -n 2 (^L x ) \ drr\d r <p(r)\ 2 (A.ll) 

Jo Jo 

where we use that the second term in ( |A.10| ) is a total derivative and that g(0) =0 due to 
the regularity condition for the metric at the origin. In addition g/r 2 — > at infinity. In our 
units (167tGjv = 1), the black hole formula is S = 4% A = 4ttA q (1 + ^f-) = 4tv(A + 8 A). 



Substituting the solution for <f>(r) for n = 1 (|A.2| ), and integrating, we get the same as in 



( |A.5| ). We checked this only numerically, but below we will show it without performing 
the explicit calculation. 
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A. 3. The two quantities are the same 

In the above computation we actually did not need to solve all the equations to the 
end in order to show that the two results are the same. 

We will rearrange the entropy formula for the scalar so that we get an expression 
that is simpler to compare with the area contribution. The lagrangian C(g^ u , (f>, V M </>) is a 
function of r. When we evaluate the gravitational action, we integrate over all coordinates 
except r. Then we first integrate over r from zero to 2tt and then multiply by n. We can 
do this both for integer or non-integer n. We denote the r integral as [logZ(n)]2,r. Then 
we have 

logZ(n)=n[logZ(n)] 27r (A.12) 



Then the entropy formula ( |1.3|) simplifies and we get 

S = - nd n {n[logZ(n)] 2lT - nlogZ(l)}| n=1 = - c\[logZ(n)] 27r | n=1 (A.13) 

And the later expression can be straightforwardly evaluated, usingi3 y/gT^ = q^ v , 



dgt- 

-2lT 



-a n [log Z matter (n)] 27r =J o drj dxdr^-g [T» u ^_ + ^d n <P + ^-^d n d^j 



»2tt 



dg^ 



dr I dxdryfgT^- 

(A.14) 
In the last line we used the equations of motion (of course |j = 0). One can check that 
the expression with the stress energy tensor gives us y^^T^ -|— | n -i = —2r] 2 rf' 2 , so 

S - S = -ifAnL x [ drrf' 2 (A.15) 



In writing QA.14J) we have only included the action of the scalar field in the computation. 
We can now show that we get the area, without using explicit expressions. This can 
be done as follows. First note that in the second line of ( |A.14|) we can use Einstein's 
equation to write T^ v in terms of the Einstein tensor, which is related to the variation of 
the gravitational action. We end up with an expression of the form 



dg» v 

drdxdryfgG^v — — 



dn 



(A.16) 



ra=l 



We define it like this because the action is logZ(n) = j AdS (^Grav — ^matter) so the field 



equations read G^ = T^. 
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This is closely related to the derivative of the gravitational part of the action. As we 
explained above we know that the gravitational part of the action has no term of order n 2 . 
Thus we know that the d n derivative of the gravitational part vanishes at order rj 2 . This 
derivative is the same as ( [A.16Q up to a total derivative term 



d n [\ogZ G ™(n)] 2 „\ n2 



= 2tt 



dxdr^fgG^y 



dg^ v 
dn 



dxy/gV^dng 



fj,r 



r=0. 



(A.17) 



The last term gives the area of the horizon, or more precisely the area of the horizon at 
order rj 2 . 

These are the same manipulations that one can do in general, but we have done all 
steps explicitly above to check that everything indeed works in situations with no U(l) 
symmetry. 

A. 4- Real scalar 

We now consider the case of a real scalar field </> = /(r)cosr, /(r) is the same as 
before but now the stress energy tensor no longer has the U(l) symmetry 



rp rpQ , rpX l2r I rp—l—i2T 



(A.18) 



And T* = T—\. The metric has the same fourier decomposition, so v and g in ( |A.8|) also 
have three fourier components. The entropy coming from the change in the area is 



b\r)2 



2tt 



drv(0) = 2ttv°(0) 



(A.19) 



where v is the constant component of v. It is easy to check that v (0) = — J drrf' 2 . 
The scalar action contributes as follows 



D|^2 



drdr(—2rf cos (fcr)) = — 2n / drrf 



f/2 



(A.20) 



So we find agreement once more, and the result is precisely half of the complex scalar. 
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Appendix B. Derivation of minimal area condition for the general case from a 
explicit calculation 

In this appendix, we obtain the minimal area condition of section 4 without using 
dimensional reduction. As in section 4, we derive this condition from requiring that the 
analytically continued solution satisfies the linearized equations of motion near r = 0. 

The metric of the n = 1 solution, which satisfies (locally) the equations of motion is 

ds 2 =dx\ + dx\ + giAdy 1 + bl t dx a )(dy j + V a dx a ) + o(r 2 ) , 

(B.l) 

9ij = h ij + x l K ij + x 2 K ij > K ~ °( r ) 

Here, yi are the directions along the surface. Now, we do the replica trick, that is, we 
change the periodicity of the r circle from 2tt to 27rn and analytically continue n to 1 + e. 
In this way, the metric will be modified to linear order in e 

ds 2 = e 2p (dr 2 + r 2 dr 2 ) + g tJ (dy l + b^dx^dy 3 + V a dx a ) + Sg (B.2) 

Where we decomposed the perturbation in a part that makes the metric smooth p = 5p = 
— elogr and a perturbation Sg that has components Sg a i, valued in all directions. For 
simplicity we work with z,z coordinates: Xi = ^-^.x^ = ^f- As a gauge condition, we 
set 5g zz = dgzz = 0. We also set 5g zz - = 0, since this variation is included in p. We require 
the perturbation, 5g a b to be periodic: Sg a b( T ) ~ ^9ab(T + 27r). 

We want to compute the linearized equation of motion 5G ZZ = 5T ZZ . In particular, 
we want to focus on the terms that can be divergent, going like 1/r near the origin. We 
find 



5R ZZ =—K z + ^(25g p z . zp - Sg. zz - V 2 5g zz ) + (regular as r ->■ 0) 



where $7 = g^Sgij and K z = K ~ 2 lK . In ( |B.3|) we neglected the terms that have yi 
derivatives because we expect them to be regular, only terms with two x a derivatives can 
contribute to this order. 

Now, since the stress energy tensor is not expected to be singular, the equations of 
motion imply that the two potentially divergent terms should cancel 



1^ = _e 

2 z ' z 

\qI^ = zi 

2 z ' z 



;b.4) 
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These are the same equations as before (|4.7|) , which are only satisfied for a periodic 
function, 5^{t) ~ 5^(t j c2tx) 1 if K z = K s = 0. Note that although the equations of motion 
are well behaved for K z = 0, the Riemann tensor diverges, as we discussed in section 2. 



This discussion is similar to the analysis in [|23| for the motion of a cosmic string. 



Appendix C. Computation of the entropy for a disk 

Here we consider a very simple example of gravitational entropy. We go through it to 
explain how one can put boundary conditions at fixed distance. 

Consider the metric ds 2 = dr 2 + r 2 dr 2 . In addition, we can have other dimensions, 
but let us assume we can ignore them. In this case, we can say that we pick an r = r c 
and we set up the boundary conditions there. We demand that the metric in the angular 
direction is 

ds 2 bdy = r 2 dr 2 (C.l) 

at the boundary r = r c . We now consider the situation with t ~ t + 2nn. We should 
consider now metrics with the same boundary condition (|C.1|) , but compatible with the 
new period. These metrics are 

ds 2 = n 2 dr 2 + r 2 dr 2 (C.2) 

We can evaluate the gravitational action for these spaces and obtain 

A 



log Z(n) 



16tvGn 



[ VdR + 2 [ K 

J Jbdy 



AG 



N 



(C.3) 



which is independent of n. Here A is the area of the transverse directions which were not 
explicitly mentioned above. Using the usual formula, we get the expected area formula for 
the entropy. 

We have included this trivial computation to explicitly show how gravity regularizes 
the divergent contribution that one normally gets in field theory. In fact, there is no 
divergence because there was no conical space in this computation!. Of course, this begs 
the question of whether the finite part of the one loop correction computed by performing a 
one loop computation around the above geometries is indeed the same as the finite part of 
the one loop corrections computed using the conical spaces that appear in the field theory 
discussion of the replica trick. 
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